Statistical description of magnetic 
domains in the Ising model 



K. Lukierska-Walasek 

Institute of Physics 
University of Zielona Gora 
ul. Z. Szafrana ^a, 65-516 Zielona Gora, Poland 

K. Topolski 

Institute of Mathematics 
Wroclaw University 
pi. Grunwaldzki 2/4, 50-384 Wroclaw, Poland 



Abstract: 

We use the Mandelbrot-Zipfs power law for the description of the inhomogenity of 
the spin system. We describe the statistical distributions of the domain's masses 
in the Ising model near the phase transition induced by the temperature. The 
statistical distribution near the critical point appears to be of the Pareto type. 

We study in this paper for the description of the phase transition the Ising 
model [1]. This model is one of the simplest models which can be used to 
describe the phase transition in ferromagnet. Besides this model is known due 
to several applications: first of all the percolation [2], then trading activity 
[3], sociophysics [4] and others. A common feature in these problems is the 
presence of two choices of the variable 5" = ±1. However all the results one 
can generalize to the model containing many values of S for instance to the 
Potts model. 

The Hamiltonian for the simplest Ising model is: 



with the sum over all neighbour pairs {z-th component) of spins. Usually it is 
assumed that the crystal lattice of ferromagnet is regular and in each site of a 
lattice the spin is locahzed with the value — 1 or — —1. Further 



Two spins i and j interact with each other by an energy —J Sf with 
—J if both spins are parallel and +J if they are oposite to each other. The 
energy needed for fliping of one spin is 2 J. 

For the simulation in this model we will use the Monte Carlo method with 
the Swendsen-Wang cluster algorithm [5]. 

In this algorithm clusters of spins are created by introducing bonds between 
neighboring spins with probability P{Sf,Sj) = 1 — exp(— ^|;), where ks is 



(1) 





J a i,j are neighbour pairs of spins 
in oposite case. 
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Boltzmann constant, AE is the energy difference needed to transform a pair 
of equal spins to a pair of opposite spins. Tlie probability is zero if spins are the 
same. All such clusters are generated and then updated by choosing a random 
new spin value for each cluster and assigning it to all spins in this cluster. 
The probability is zero if spins are the same. Starting with the simmulation 
having the random distribution of half of the spins up and half down and using 
Swendsen-Wang algorithm with low temperature one sees growing domains, in 
which spins are parallel. We have two kind of domains: with spins up and 



there appears an infinite domain in the limit L ^ oo where L denotes linear 
size of the system with one of the spin states being chosen. 

The phase transition appears in the critical point T = Tc. The difference 
M between a number of spins up and down is proportional to the magnetisa- 
tion and near critical point vanishes as {T — Tc)^, where for dimension d — 2, 
(3—1- The correlation length ^ ~ IT—Td". The magnetisation is proportional 
to ^"^^^ . In a finite system in critical temperature Tc one can replace ^ by L, 
hence M ~ L'^-PI'' = with the fractal dimension D ^d- /S/u {d<4). 

The simulation data were collected on square lattices of linear size L = 
500 and 1000. A total of 10000 Monte Carlo (MC) time steps were used for 
equilibration. The value of MC time steps required for equilibration have been 
estimated from the energy time series, which is a common practice for cluster 
algorithms. 

The main goal of this paper is the statistical description of the simple mag- 
netic system when we approach the critical point of phase transition induced 
by the temperature. 

For this purpose we will consider Mandelbrot-Zipf's power law [6]: 



In our case x is the number of spins up or down in the domain (domain mass) , 
k denotes the rank order of the domain mass x. (The greatest cluster has rank 
1, smaller rank 2 and so on.) 



with spins down. At last at the temperature Curie T = T^ 




X — k 
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Bouchaud [5] pointed out the strong correlation between Mandelbrot-Zipf 's 
power law and the inhomogenity of the system: the slope (in log-log) of the 
straight line is determined by — ^ and characterize the inhomogenity of the 
physical structure of the system (— ^ = tga, where a the angle of the slope). 
The inhomogenity of the system means that its structures become fractal and 
more hierarchical. 

In our considerations we shall concentrate on the sequence of random vari- 
ables X called here the [/x]- variables. These variables are distributed according 
to the distribution of the appearance of a cluster with mass and probabil- 

-j.n 

ity p[x) which dacays as ^ttti where xq is the typical scale. The index 
appearing in the tail of distribution p[x) is a critical exponent. The main 
property of [p\ variable is that all its moments ruq =< x'^ > with q > /i are 
infinite. 

We are going to connect the statistics of domain masses with the pro- 
cess of approaching the critical point. When we start to advance from the 
paramagnetic phase to the critical point {T — > Tc) as we see on Fig (1) 

Insert Fig (1) 

Fig 1. Log- log distribution of the domain's masses x versus the rank order 
index A; (L = 1000). 

1/T = 0.25: ln(ln(k)) = 4,01287428 - 0,148950538*ln(k) 
1/T = 0.35: ln(ln(k)) = 5,55609713 - 0,21268554*ln(k) 
1/T = 0.4: ln(ln(k)) = 6,89830264 - 0,261512561*ln(k) 

1/T = 0.44068: ln(ln(k)) = 10,4914484 - l,014818*ln(k) 

The angle between strait hues representing (in log-log) the Mandelbrot-Zipf 's 
inverse power law and the rank axis increase when // > 1 and T — > Tc- 

For that case we observe the growing domains, their structure become 
more fractal (loss of an oval) and more hierarchical - the inhomogenity of the 
system increase. The distribution of magnetisation of a whole system has the 
usual Gauss form. At high temperature correlations between spins are short 
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ranged in the whole high temperature region M ~ = L for d = 2, see 
Fig. 2. 

Insert Fig (2) 

Fig 2. Test M ~ L for T = 5 T^. 

At last at T = Tc and fJ, ^ I, when one domain covers the whole lattice and the 
other spin orientation is restricted to the small clusters or isolated single spins 
within the domain. In this case the distribution of the total magnetisation is 
centered at mL'^ or —mL'^, where m is the remnant magnetisation and it 
also has Gaussian form, as it is well known. 

In the case when // ?a 1 {T — Tc) from fig (1) we see the strait hne 
representing (when a — ac, tg etc = ^ ~ 1) the Mandelbrot law describing 
the phase transition. The highest point denotes in the critical temperature the 
domain which covers the whole lattice. 

Fig (3) represents the histogram of the domain masses at critical point - the 
probability of the appearence of the cluster with mass x. 

Insert Fig (3) 

Fig 3. Histogram of the domain's masses at critical point T — Tc {L — 500) 

The distribution p[x) of domain mass is like has Pareto tail ^1+7:, where xo 
denotes a typical scale, with /i 1. 

When II > 1 distribution p{x) is without power-law tail, which becomes 
truncated, see Fig 4. 

Insert Fig (4) 

Fig 4. The histogram of the domain's masses (/3 = ^ = 0, 3; L — 500) 

This results are in the agreement with standard percolation theory [2] and the 
paper of Janke and Schakel [6] because the distribution of domains with the 
mass X takes a general form 

Qx ~ x~'^ exp{—Qx) 
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where r is the entropy factor (r = + 1), - (T - TJ^, if T ^ T^. When 
T = Tc, we get © = and we have ~ x~'^. 

The statistical description of discontinuous metal films on dielectric sub- 
strates were analysed on the basis of experiments by Dobierzewska-Mozrzymas 
at all [8] and the distributions of local fields intensities in metal dielectric sys- 
tem was investigated by Liberman et al. [9]. Unfortunately the analogous 
experiments on ferromagnetics are not known. 

The main result of our paper is to show the conection beetwen the 
Mandelbrot-Zipf 's law and the statistics of the domain masses in the Ising 
model. The domain masses in the Ising model fulfil the Mandelbrot- Zipf's 
inverse power law and when we approach the phase transition in this model 
the distribution of the domain masses appears to have the Pareto tail. This 
model is a such one which represents the system in which in the critical point 
the lenght scale diverges and leaves the system in self similar state. That 
feature denotes a fined-tuned criticality which should be contrasted with self- 
organized criticality. In such a case the system spontaneously evolves towards 
scale-invariant states and one raise the problem raises problem of the univer- 
sality of the renormalized coupling constant at critical point. This conclusion 
is in agreement with the paper of Hilfer [10]. 
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